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We present a general scheme to determine the loss-free adiabatic eigensolutions (dark-state po- 
laritons) of the interaction of multiple probe laser beams with a coherently driven atomic ensem- 
ble under conditions of electromagnetically induced transparency. To this end we generalize the 
Morris-Shore transformation to linearized Heisenberg-Langevin equations describing the coupled 
light-matter system in the weak excitation limit. For the simple lambda-type coupling scheme the 
generalized Morris-Shore transformation reproduces the dark-state polariton solutions of slow light. 
Here we treat a closed-loop dual-!/ scheme wherein two counter-propagating control fields generate 
a quasi stationary pattern of two counter-propagating probe fields - so-called stationary light. We 
show that contrary to previous predictions, there exists a single unique dark-state polariton; it obeys 
a simple propagation equation. 



I. INTRODUCTION 

The adiabatic interaction of probe light pulses with 
coherently driven three-level systems under conditions of 
electromagnetically induced transparency (EIT) [ij can 
most convenientlybe described in terms of so-called dark- 
state polaritons [4]. These are the eigensolutions of the 
Raman interaction of probe and control fields with the 
atomic ensemble in the absence of atomic losses. They 
provide a simple theoretical framework for phenomena 
such as slow light Q , the storage and retrieval of coherent 
light pulses in atomic ensemble [3, Q and quantum mem- 
ories for photons 0, Q . Most importantly they also fully 
incorporate the essential physics of adiabatic pulse prop- 
agation for time-varying control fields together with the 
associated coherent transfer of excitation and quantum 
state from light to matter and vice versa. This transfer 
cannot be understood in terms of electromagentic field 
equations alone. 

Recently it has been shown that the simultaneous pres- 
ence of two counter-propagating control fields of compa- 
rable strength can lead to a quasi-stationary pattern of 
slow light consisting of two counter- propagating probe 
field components 0, la, HI El ■ Such stationary light pulses 
hold particular interest as examples of efficient nonlinear 
optical processes. In contrast to the three-level coupling 
scheme, there exists no unique description of stationary 
light in terms of dark-state polaritons. Reference H, [llj] 
introduced separate polariton fields for forward and back- 
ward propagating modes. However, as we will show here, 
these are not the quasi loss-less eigensolutions of the sys- 
tem. 

We here determine the adiabatic eigensolutions of a 
coupled multi-level system using a generalization of the 
stationary-light scheme. As with simpler systems, the 
eigensolutions are superpositions of collective atomic op- 
erators and operators of the electromagnetic field. Our 
procedure relies on the possibility of casting the lin- 



earized Maxwell-Bloch equations into a matrix differen- 
tial equation similar to that used in treating the time- 
dependent Schrodinger equation (TDSE) describing co- 
herent excitation of a multi-level atom. For that equation 
there exists a well-known formalism to determine possi- 
ble non-decaying eigenstates, the Morris-Shore transfor- 
mation fl2| . Generalizing this transformation to the lin- 
ear Heisenberg-Langevin equations of the coupled light- 
matter system we can easily determine their loss-free 
eigensolutions, when they exist. We apply this proce- 
dure to the stationary-light equations to find a dark-state 
polariton solution, and to derive its equation of propaga- 
tion. 



II. THE DUAL-V STATIONARY-LIGHT 
SCHEME 

Let us consider one-dimensional propagation in the z 
direction of two counterpropagating probe fields, circu- 
larly polarized orthogonally, having electric fields opera- 
tors E±(z,t) and propagation vectors k±. These inter- 
act with an ensemble of stationary four-state atoms via 
linkages between the ground state \g) and excited states 
\e±). A second pair of counterpropagating fields, also cir- 
cularly polarized orthogonally and termed control fields, 
link these excited states to a fourth state \s), through 
interactions characterized by Rabi frequencies Q±(z,t) . 
Figure [T] illustrates the linkage pattern of the four fields 
with a four-state atom. The control fields are assumed to 
be sufficiently strong that they are not influenced by the 
atoms; unlike the probe fields, they are not considered as 
dynamical variables. 

We assume that the control fields generate 
electromagnetically-induced transparency (EIT) for 
the probe fields, by requiring that their Rabi-frequencies 
are locked to those of the control fields. Thus if 
= the control fields form a standing- wave 
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FIG. 1: The dual-V linkage pattern shown together with 
the depicted field configuration of counter-propagating con- 
trol (n±) and induced probe fields (E±) leads to stationary 
probe pulses. 



pattern and the group velocity of the probe fields vanish; 
these form a standing-wave pattern known as stationary 
light 0,0. 

A stationary-light probe-field pulse can be created 
from a propagating-light probe field by adiabatically 
changing the relative intensity of the two control fields 
from zero to unity, thereby storing the probe pulse in 
a collective spin excitation [To| . The stored excitation 
can then be retrieved as a stationary pulse by using two 
counter-propagating control fields of equal intensity. 

Stationary-light pulses were suggested and experimen- 
tally demonstrated for the first time for a A-type three- 
state configuration by Lukin and co-workers @, Q. In 
contrast to the system discussed here, there the counter- 
propagating fields, whether probe or control field, acted 
upon the same transitions. A simple theoretical descrip- 
tion of the A-type system of [9( is possible, by neglecting 
components of the atomic spin coherence that are rapidly 
oscillating in space. Such a secular approximation is only 
justifiable for a vapor at sufficiently high temperature, 
such that thermal motion of the atoms leads to a rapid 
dephasing of the rapidly oscillating components 0, [Tl| . 
Here we avoid this dephasing approximation and consider 
the 4-level generalization shown in Fig. [T] This model 
also produces stationary light but does not require any 
secular approximation. 

We assume that the control fields have carrier fre- 
quencies Wjf and wave-vectors fc^f = ±fc( c ) directed 
along the z axis. We take the carrier frequencies and 
wave-vectors of the probe-field components to be uj^ 



and k± = ±k( p \ again along the z axis. We replace 
the basic parameters of the control fields £l±(z, t) and 
the operators E±(z,t) of the probe fields by quanti- 



ties n±(z,t) = ft±(z,t)e- i(fc ± 

E±(z,t)e~ i[k ± )z ~ u '± )t) that vary slowly with z and t. 
We further require four-photon resonance and phase- 



J ± 4) and £±(z,t) 



matching, as expressed by the constraints 
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Electromagnetically induced transparency requires that 
the control and probe fields for both the forward and 
the backward propagation direction are in two-photon 
resonance with the atomic system, i.e. 
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where tu sg is the transition frequency between the ground 
state \g) and the spin state |s). 

To describe the atomic system we use collective atomic 
variables ^y, v {z,t), generalizations of the traditional 
single-atom transition operators = \n){v\ between 
states \v) and | /j.) 0, 0]. We take the two probe fields 
to have the same field-atom coupling strength, denoted 
g, and we denote by N the one-dimensional density of 
atoms, taken to be uniform. Then the interaction with 
the light fields is described, in the rotating-wave approx- 
imation, by the interaction Hamiltonian 



H int = hJ2 / '!•-(-. g &, - fij 



i(fe (c) z-w' c) t) ~ 



(3) 



The rapidly oscillating exponents in §3§ can be eliminated 
by introducing slowly varying amplitudes of the atomic 
variables 



,(z,t) e 
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T e±g (z,t) = o- e±g (z,t) e n " 
a sg (z,t) = a sg {z,t)e 



(4) 
(5) 
(6) 



In contrast to the original stationary- light scheme 0, [T(| , 
the conditions ((T|) allow the introduction of a slowly- 
varying amplitude of the ground state coherence a gs 
without a secular approximation. The coherent inter- 
action Hamiltonian can now be rewritten as 



Here A± 



,(p) 



+gVN£ j a ejg + h. a. . (7) 

are the one-photon-detunings 
of the forward (+) and backward (— ) propagating probe 
fields respectively. 

The atomic system is also subject to losses and deco- 
herence. These we treat by a set of Heisenberg-Langevin 
equations for the collective atomic operators and the op- 
erators for the probe fields [l5[ ■ We assume that the con- 
trol fields are sufficiently strong to remain undepleted - 
they are not dynamical variables - and we consider only 
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the linear response of the atoms to the probe fields. In 
this limit we treat the probe fields as perturbations in 
the Heisenberg-Langevin equations. 

In zeroth order all atoms will be optically pumped into 
the ground state \g). The only nonzero atomic variable 

in this limit is a g °J = 1; other operators vanish, affl = 0. 

In first order one finds a linear set of Heisenberg- 
Langevin equations for the atomic variables. These sep- 
arate into two uncoupled sets. Those relevant for this 
work are 



<9 t <7 



tO ge± 



dt&gs 



[ lA± +7 S e±] o>± 

igV~N£± + iQ±a gs - 
i£l* + a ge+ + ifl*_a ge _, 



F 



(8) 
(9) 



where the quantities "f ge ± parametrize the losses and de- 
coherence rates and the F ge ± are Langevin noise forces 
EE 

UH associated with the decays. In eq. © we have 
assumed that the ground-state coherence a gs is stable. 
This is appropriate for hyperfine-states in cold and hot 
atomic vapours for which the coherence lifetimes, mainly 
caused by elastic collisions, are of the order of millisec- 
onds HE H- As a consequence of this stability the 
equation includes no corresponding Langevin noise oper- 
ator. 

The Langevin noise operators FA(t) are necessary to 
preserve the commutation relations of quantum variables. 
We assume, as is customary, that the associated decay is 
exponential in time, so that the noise operators are delta- 
correlated in time, 



(F A (t)F B (t')) = D AB (t)5(t-t r ). 



(10) 



The diffusion coefficient D AB (t) appearing here is re- 
lated to the dissipative part of the dynamics through the 
dissipation-fluctuation theorem [l5|, fill ] 



D AB (t) 



w 'A(t)B(t) 
at \ / loss 



B(t) 



loss 



-(A(t)-B(t) 



loss 



One easily verifies that, because the quantum noise orig- 
inates from spontaneous emission processes, the relevant 
diffusion coefficients of the present system are propor- 
tional to the excited-state populations. However, in the 
linear response limit considered here this population is 
negligible and thus we can ignore the Langevin noise 
terms altogether. 

To complete the description of the atom-field system 
we require equations for the probe fields. These we take 
to be wave equations for the slowly-varying probe field 
amplitudes 



(11) 



[dt ± cd z ] £ ± = igVNa : 



These, together with the atomic equations © and © 
form the set of self-consistent set of Maxwell-Bloch equa- 
tions which are the basis of the following considerations. 



We note that equations ([8]) and ([9]) are formally iden- 
tical to the time-dependent Schrodinger equation (or 
TDSE) for the A-system after the secular approximation 
(see e.g. eq. (9) and (10) in [l(|). No such approximation 
has been used here, however. Furthermore, in contrast 
to the system studied in 0, [l(| the single-photon detun- 
ings A± can be chosen to be different for the forward 
and backward direction, which adds another degree of 
control. 

We now assume that the slowly varying amplitudes 
of the control field Rabi-frequencies are constant in 
space, as is the atom density N. We introduce a spa- 
tial Fourier transformation, through variables X(z,t) 
I dk exp(-ikz) X(k,t), thereby reducing the Heisenberg 
equations and the shortened wave equation to coupled 
differential equations in time alone. We write these in 
matrix form as 



dt 



X = -MX 



(12) 



where the column vector of dynamical variables has the 
elements X T = {£+, a gs , a ge+ , o ge _ }, and the matrix 
of (slowly varying) coefficients is 



(13) 



To simplify typography we have introduced the symbol 

(14) 
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The linkage pattern of these equations is shown in Fig. 
[2] It corresponds to the M-type linkage found with the 
TDSE for light fields interacting with atoms, as was ex- 
tensivly studied in the literature In the following 

we draw on this analogy with the conventional TDSE 
to identify the loss-free adiabatic eigensolutions of the 
system through use of the Morris-Shore transformation, 
a method widely used to find adiabatic dark states of 
the interaction of atomic multi-level systems with near- 
resonant laser fields. 



III. THE MORRIS-SHORE TRANSFORMATION 

Let us first summarize the basic properties of the 
Morris-Shore transformation [ijj , before considering the 
two examples most relevant for us. 

Consider a first-order ODE, of the form (fT2")) , in which 
X is an iV-dimensional column vector that separates into 
two classes of variables: a set A of Na variables and a set 
B of Nb variables. We suppose that the N x N matrix H 
has elements only between the A and B sets, not within 
them. X and H therefore have the forms 



X = 



x A 
x B 



H 



V 

V+ 



(15) 
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where V is an N A x Nb matrix. We further assume 
that initially (at t — 0) all variables of the B set van- 
ish, Xb(£ = 0) = 0. Such equations can, by means of 
the Morris-Shore (MS) transformation Y = M T XM fl^ . 
be rewritten as a set of independent equations involv- 
ing only pairs of variables (one being a combination of 
the A set variables, the other a combination from the B 
set) or unlinked single variables. The number of unlinked 
("dark") variables N D is the difference N D = \N A -N B \- 
In the new basis the equation of motion reads 



dt 



Y = -iH 



(16) 



where the matrix of coefficients has the structure 



•MS 
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(17) 



Here the are matrices of dimension 2x2 



,,(i) 





yU) 



.j(J) 





(18) 



There are mm(N A . 
"bright" variables, 
main constant, i.e. 



Nb) of these, each linking pairs of 
The remaining "dark" variables re- 



dt 



Y n = 0. 



(19) 



For the A and dual-V^ linkages considered here the A 
set includes variables representing both field and atom 
degrees of freedom. The resulting MS dark variables, 
combining field and atom properties, are dark-state po- 
laritons. 



A. The A-system 

For the A-system (N = 3 with N A = 2,N B = 1 and 
hence Nd — 1 dark state) the original matrix, when ex- 
pressed in a basis {Xi, X 2 , X 3 } , has the form (apart 
from reordering) 



Vi 
V 2 

Vi v 2 



(20) 



The single dark variable, constructed entirely from A 
variables, is 



Y D = [V 2 X l - V X X 2 ] IH 



(21) 



where M is a normalization factor. Note that Yd remains 
a "dark" variable even when a finite diagonal element is 
added to the third row of the matrix in (f2"0"l) . 



ge- 



ge+ 



1/2 ^ ^ 1/2 



a, 



kc 
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■ £_ 



gs 



-kc 



FIG. 2: The M-scheme corresponding to the dual-!/ coupling 
scheme depicted in Fig. [T] The matter and field amplitudes 
are treated as basis states of the new system. 



To connect this formalism with that of dark-state po- 
laritons we let X represent the amplitudes of the wavevec- 
tor of the atomic A system, with X\ and X 2 being the 
amplitudes of the stable (i.e. non-decaying) lower levels. 
Then Yd is the well-known dark state of that system. 



B. The M-system 

The M-system (N = 5 with N A = 3, = 2 ) has also 
only one dark variable. The original matrix has the form 
(apart from reordering) 



Vi 

v 2 

v 3 14 

Vi v 3 

v 2 v 4 



The single dark variable is 

Yd = [V 2 V 3 X 1 + VxViXz - V^X^ /M 



(22) 



(23) 



where Af is a normalization constant. As with the A- 
system, Yd remains a "dark" state if diagonal elements 
are added to the two lowest rows of the matrix. 



IV. APPLICATION OF THE MORRIS-SHORE 
TRANSFORMATION TO THE DUAL-V SYSTEM 

We now apply this mathematical tool to actual phys- 
ical problems. The dual-V system depicted in Fig [TJ can 
be described in the Heisenberg picture as the M-system 
shown in Fig. [2] Applying the Morris-Shore transfor- 
mation one can immediatly construct the dark variable. 
Neglecting the (±fec) diagonal elements we find 



Yd = n+gVN£ + + n_gVN£ 



g 2 Na g 



IN. (24) 
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To express this more simply we define mixing angles by 
means of the relationships tan 2 8 — g 2 N/il 2 and tan 2 ip = 
^ 2 _/^+, where fl 2 = + fi 2 . The construction then 
reads 



Yd = ^cos ip£+ + sin ip£- j cos 8 



a as sm ( 



(25) 



This is the Fourier transform of the dark-state polariton 
(DSP) for a dual-y-system. To obtain the adiabatic limit 
we did set kc — > in eqs. (fTB"]) , thereby neglecting all 
spatial variations; that limit produces the equation 



d_ 

dt 



Yd 



0. 



(26) 



In addition to the usual adiabatic condition £l c ffT 1, 
where fL c g — y / g 2 N + O 2 and T characterizes the pulse 
duration, one must impose a further adiabatic condition 
on the spatial Fourier frequencies to justify the negelct of 
the terms ±kc. As has been discussed in detail in Ref. 6] 
the corresponding condition is 



Lp ^S> \J -Z^abs \ I ~ > 

V 7 



(27) 



where i a b s = c-f/g 2 N is the absorption length of the 
medium, and L p is the characteristic length scale over 
which changes occur in the amplitude of the dark-state 
polariton in the medium. 

The adiabatic limit, eq. (|26D . does not account prop- 
erly for propagation effects [13|, e.g. it neglects conse- 
quences of nonzero group velocity. In order to describe 
the propagation dynamics of the DSPs one needs at least 
the lowest order corrections to the spatial Fourier compo- 
nents. These can be calculated in a very simple pertur- 
bative way. The introduction of nonzero kc ^ amounts 
to adding the terms {— kc, +kc, 0, 0, 0} to the diagonal 
elements of H, leading to the equation 



dt 



-kc 0" 

+kc 









MY. (28) 



From this we deduce the dispersion relation for the dark- 
state polariton 



■(&) = kC\ +k 2 C 2 . 



(29) 



The coefficients appearing here are obtained as the first 
two non-vanishing corrections to the eigenvalue of the 
dark-state polariton due to the terms proportional to 
±kc. 



C\ = v gr COS 2lfi 

C2 = VgrL^bs ^ i ) (siir 2y — cos" 2 r -siii 1 1 



(30a) 



(30b) 



Here Vg r — c cos 2 8 is the group velocity of standard slow 



light [2| and we have assumed that T ge+ = T ge _ = T = 
7 + iA. We see that the effective velocity of the dark- 
state polariton in the stationary light case is then given 
by 



dto 2 n 

v = — = ccos 9 cos 2ip 
dk 



(31) 



and that the DSP gains an effective complex-valued mass 

(32) 

— - i \ (sin 2 2tp + cos 2 2ip sin 4 8) . 



J_ 1 d 2 UJ 

m* h dk 2 

_ A-K Vg r L abs 



m v r 



A 



Here A = 2Tr/k^ is the wavelength of the probe field, 
m is the atomic mass, and mv rec — Hk^' the recoil mo- 
mentum of a probe-field photon. In general the effective 
mass of the DSPs of stationary light is much smaller than 
the real mass of the atoms, even if the group velocity v gr 
is of the order of the recoil velocity v rec , because the 
wavelength A is much smaller than the absorption length 

^abs- 

From the above considerations we obtain a very simple 
equation of motion for the dark-state polaritons including 
the lowest order corrections. The equation for Fourier 
components reads 



d 
dt 



ikv 



2m* 



0. 



Transformed to real space this equation becomes 



d_ 
dt 



d 

oz 



h d 2 



2m* dz 2 



= 0. 



(33) 



(34) 



The effective group velocity of the DSP, eq. (|3ip . is fully 
determind by the mixing angles. The DSP comes to rest, 
V = 0, when at least one of the two trigonometric func- 
tions, cos# or 008 2^, vanishes. It is interesting to note 
that the sign of v can be changed by altering the mixing 
angle tp, i.e. by changing the ratio of the Rabi frequen- 
cies of the two counter-propagating control fields. If both 
these fields have the same strength then v — and the 
field pattern has no net motion. If the forward (back- 
ward) propagating control field is stronger than the other 
one, the net motion of the probe field pattern is also in 
the forward (backward) direction. 

One also notices that when cos 2^ = the first-order 
term C\ vanishes, and thus the dynamics is dominated 
by either a Schrodinger-type evolution or a weak diffusive 
spreading of the probe field, depending on the relative 
size of imaginary and real parts of m* . On the other 
hand, if only one of the two control lasers is present then 
cos2(/3 = ±1 and we recover the equation of motion for 
the slow light polaritons found in 0]. 

It should also be noted that as soon as one knows the 
form of the dark state one is able to calculate all other 
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states of this system, also termed bright-state polaritons. 
This can be done by defining mutually orthogonal eigen- 
solutions, starting with the dark state. One is able also 
to derive all non-adiabatic couplings arising e.g. from 
time dependent mixing angles, leading to first order cor- 
rections to the adiabatic approximation. 

Let us finally address some limiting factors of the con- 
sidered system/approach in particular phase fluctuations 
of the control fields, collisional dephasing of the ground- 
state transition and Doppler broadening due to atomic 
motion. In the above discussion we have disregarded pro- 
cesses which destroy either the amplitude or the phase of 
the ground-state coherence a gs . While the amplitude of 
ground-state superpositions is usually quite robust, its 
phase can be destroyed by two processes: atomic col- 
lisions and fast phase fluctuations of the coupling laser 
[ijj ■ If the laser phase fluctuations are sufficiently slow, 
such that adiabaticity holds true, probe and control fields 
are phase correlated and their difference phase does not 
change Phase destroying collisions as well as fast 

laser phase fluctuations lead to a decay of dark-state po- 
laritons into bright polaritons which are subsequently ab- 
sorbed. Thus the concept of the dark polaritons is only 
useful if the latter process is sufficiently slow on the time 



scales of interest. 

Let us secondly introduce the Doppler width of 
the \g) — | e±) -transitions and the factor r± = (w e ±g — 
uj e ± s )/uj e ±g. Using these we note, that if the follow- 
ing conditions: Q 2 » -f gs A£>, Q 2 » Aor± and 
£1 2 >> (7 ge ± + lgs±)A.Dr± are met, the influence of 
Doppler broadening onto the absorptive and dispersive 
properties of a A-type three-level system is negligible as 
long as the Doppler-free configuration of co-propagating 
control and probe fields is used [n| ■ The dual-V station- 
ary light scheme we have introduced here is essentially a 
double A-scheme with the two linkages coupling to the 
same coherence a gs . Both A-type sub-systems work in 
the necessary Doppler-free configuration. If the lower 
states \g) and \s) are hyperfine levels, the above con- 
ditions are practically always fulfilled. Hence, Doppler 
broadening is not a limiting factor here. 
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